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Abstract. Previous work has related the Galactic bar to structure in the local stellar ve- 
locity distribution. Here we show that the bar also influences the spatial gradients of the 
velocity vector via the Oort constants. By numerical integration of test-particles we sim- 
ulate measurements of the Oort C value in a gravitational potential including the Galactic 
bar. We account for the observed trend that C is increasingly negative for stars with higher 
velocity dispersion. By comparing measurements of C with our simulations we improve on 
previous models of the bar, estimating that the bar pattern speed is Ob/f2o = 1-87 ± 0.04, 
where Q. is the local circular frequency, and the bar angle lies within 20° < (/>o < 45°. We 
find that the Galactic bar affects measurements of the Oort constants A and B less than ~ 2 
km s _I kpc~' for the hot stars. 



1. Introduction 

The Galaxy is often modeled as an ax- 

isymmetric disk. With the ever increasing 

proper motion and radial velocity data, it is 

now apparent that non-axisymmetric effects 

can not be neglected (nonz ero Oort constant 

C, lOlling & Defined I2003I hereafter OD03 : 

i — ii 1 

nonzero vertex deviation, Famaev et al. 2005; 

asymmetri es in the loca l velocity distribution 
of stars, IDehned [l998h . It is still not clear, 
however, what the exact nature of the perturb- 
ing agent(s) in the Solar neighborhood (SN) 
is(are). Possible candidates are spiral density 
waves, a central bar, and a triaxial halo. 

Due to our position in the Galactic disk, 
the properties of the Milky Way bar are hard 
to observe directly. Hence its parameters, such 
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as orientation and pattern speed, have only 
been inferred indirectl y from observations of 
the inner Galax y (e.g., [Blitz & Spergell lT99lt 
IWeinberdll992h . 

However, the bar has also been found to 
affect the local velocity distribution of stars. 
HIPPARCOS data revealed more clearly a 
stream of old disk stars with an asymmetric 
drift of about 45 km s _1 and a radial velocity 
u < 0, with u and v positive toward the Galactic 
center and in the direction of Galactic rota- 
tion, respectively. This agglomeration of stars 
has been dubbed the 'Hercules' stream or the 
'M-anomaly' . The numerical work of IDehned 
d 19991 120001) and |Fu! d200ll) has shown that 
this stream can be explained as the effect of the 
Milky Way bar if the Sun is placed just out- 
side the 2: 1 Outer Lindblad Resonance (OLR). 
Using high-re solution spectra of n earby F and 
G dwarf stars, Bensbv et al. (2007) have inves- 
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tigated the detailed abundance and age struc- 
ture of the Hercules stream. Since this stream is 
composed of stars of different ages and metal- 
licities pertinent to both thin and thick disks, 
they concluded that a dynamical effect, such as 
the influence of the bar, is a more likely expla- 
nation than a dispersed cluster. 

Assuming the Galactic bar affects the shape 
of the distribution function of the old stellar 
population in the SN, an additional constraint 
on the bar can be provided by considering the 
values derived for the Oort constant C. In other 
words, in addition to relating the dynamical in- 
fluence of the bar to the local velocity field, C 
provides a link to the gradients of the velocities 
as well. The study of OD03 not only measured 
a non-zero C, implying the presence of non- 
circular motion in the SN, but also found that C 
is more negative for older and redder stars with 
a larger velocity dispersion. This is surpris- 
ing as a hotter stellar population is expected to 
have averaged properties more nearly axisym- 
metric and hence, a reduc ed value of \C\ (e.g., 
iMinchev & Ouillenll2007l; hereafter Paper I). 

It is the aim of this work to show it is 
indeed possible to explain the observationally 
deduced trend for C (OD03) by modeling the 
Milky Way as an exponential disk perturbed 
by a central bar. By performing this exercise, 
we provide additional constraints on the bar's 
pattern speed. Models for the structure in the 
bulge of our galaxy are difficult to constrain 
because of the large numbers of degrees of 
freedom in bar models. Thus future studies of 
structure in the Galactic Center will benefit 
from tighter constraints on the parameters de- 
scribing the bar, such as its pattern speed and 
angle with respect to the Sun. 



2. The Oort constants 

We can linearize the local velocity field (e.g., 
Paper I) about the LSR and write the mean ra- 
dial velocity and longitudinal proper motion 
Ji l as functions of the Galactic longitude / as 

V 4 = K + Asm(2l) + Ccos(2l) (1) 
d 

Ji, = B + A cos(20 - C sin(2/) 



where d is the average heliocentric distance of 
stars, A and B are the usual Oort constants, and 
C and K are given by 

2Cs -»- + f- l - d 4± (2) 

r Or r o<f> 

2K = +- + 1 r + -^r- (3) 
r Or r oq> 

Here r and <p are the usual polar coordinates 
and = vo + v, where Vo is the circular veloc- 
ity at the Solar radius, ro. In this work we pri- 
marily consider a flat rotation curve (RC; see, 
however $5]), hence the derivatives of in the 
above equations are identical to the derivatives 
of v. C describes the radial shear of the velocity 
field and K its divergence. For an axisymmet- 
ric Galaxy we expect vanishing values for both 
C and A0- Whereas C could be derived from 
both radial velocities and proper motions, K 
can only be measured from radial velocities, in 
which case accurate distances are also needed. 

A problem with using proper motions data 
has been described by OD03. The authors 
present an effect which arises from the longi- 
tudinal variations of the mean stellar parallax 
caused by intrinsic density inhomogeneities. 
Together with the reflex of the solar motion 
these variations create contributions to the lon- 
gitudinal proper motions which are indistin- 
guishable from the Oort constants at < 20% 
of their amplitude. OD03 corrected for the 
'mode-mixing' effect described above, using 
the latitudinal proper motions. The resulting C 
is found to vary approximately linearly with 
both color and asymmetric drift (and thus mean 
age) from C ~ km s _1 kpc 1 for blue sam- 
ples to C » -10 km s _1 kpc 1 for late-type 
stars (see Figs. 6 and 9 in OD03). Since C is re- 
lated to the radial streaming of stars, we expect 
non-axisymmetric structure to mainly affect 
the low-dispersion population which would re- 
sult in the opposite behavior for C. In Paper 
I we showed that spiral structure failed to ex- 
plain the observed trend of C. 

Note that the Oort constants are not con- 
stant unless they are measured in the SN. Due 

1 Note, however, that C and K would also be zero 
in the presence of non-axisymmetric structure if the 
Sun happened to be located on a symmetry axis. 
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Parameter 


Symbol 


Value 


Solar neighborhood radius 


'o 


1 


Circular velocity at ro 


Vo 


1 


Radial velocity dispersion 




0.05v or0.18v 


cr u scale length 


r<r 


0.9r 


Disk scale length 


r p 


0.37r 


Bar strength 


f% 


-0.012 


Bar size 


I'b 


0.8r cr 



to non-axisymmetries they may vary with the 
position in the Galaxy. Thus, the Oort con- 
stants have often been called the Oort func- 
tions. 

3. The simulations 

We perform 2D test-particle simulations of 
an initially axisymmetric exponential galactic 
disk. To reproduce the observed kinematics of 
the Milky Way, we use disk parameters con- 
sistent with observations (see Table [TJ). A de- 
tailed description of our disk model and sim- 
ulation technique can be found in Paper I. We 
are interested in the variation of C with color 
(B - V), and asymmetric drift v a . We simulate 
variations with color by assuming that the ve- 
locity dispersion increases from blue to red star 
samples. For a cold disk we start with an initial 
radial velocity dispersion <x„ = 0.05vo whereas 
for a hot disk we use cr u = 0.18vo. The back- 
ground axisymmetric potential due to the disk 
and halo has the form <t>o(r) = log(r), corre- 
sponding to a flat RC. 

3. 1 . The bar potential 

We model the nonaxisymmetric potential per- 
turbation due to the Galactic bar as a pure 
quadrupole 

<D b = A b (e b ) cos[2(0 - O b /)] x (4) 

Here A b (e b ) is the bar's gravitational poten- 
tial amplitude, identical to the same name pa- 



rameter used by iDehnenl d2000l) : the strength 
is specified by e b = —a from the same pa- 
per. The bar length is r b = 0.8r cr with r cr 
the bar corotation radius. The pattern speed, 
Q b is kept constant. We grow the bar by lin- 
early varying its amplitude, e b , from zero to 
its maximum value in four bar rotation peri- 
ods. We present our results by changing Q b and 
keeping ro fixed. The 2:1 outer Lindblad res- 
onance (OLR) with the bar is achieved when 
Q b /Qo = 1 +k/2 * 1 .7, where k is the epicyclic 
frequency. We examine a region of parame- 
ter space for a range of pattern speeds placing 
the SN just outside the OLR. For a given pat- 
tern speed one could obtain the ratio ro/roLR 
through ro/roLR = ph/Qou t ~ Oh /1-7. 

In contrast to IDehnenl fe OQO) and simi- 
lar to iFuxl d200ll) and iMuhlbauer & Dehnenl 
(120031) . we integrate forward in time. To al- 
low for phase mixing to complete after the bar 
reaches its maximum strength, we wait for 10 
bar rotations before we start recording the po- 
sition and velocity vectors. In order to improve 
statistics, positions and velocities are time av- 
eraged for 10 bar periods. After utilizing the 
two-fold symmetry of our galaxy we end up 
with ~ 5 x 10 5 particles in a given simulated so- 
lar neighborhood with maximum radius af max = 
ro/10. 

4. Results 

First we show that we can reproduce the results 
of Dehnen and Fux for the Hercules stream. 
We present a simulated u — v velocity distri- 
bution for Q b = 1.87Qo, (po = 35°, and a sam- 
ple depth d = r /40 in Fig. [TJ The initial ve- 



192 



Minchev, Nordhaus, and Quillen: Milky Way bar pattern speed 




0,4 -0.2 0.0 0.2 0.4 
u/v > GC 



Fig. 1. Simulated u - v distribution for a bar pattern 
speed Q. b = 1.87f2 , bar angle O = 35° and a sam- 
ple depth d = ro/40. The initial velocity dispersion 
is cr„ = 0. 18v and the bar strength is e b = -0.012. 
Contour levels are equally spaced. The clump iden- 
tified with the Hercules stream is clearly discernible 
in the lower left portion of the plot. This figure is in 
agreement with the observed fo(u, v) (Dehnen 1998; 
Fuxll200l | ) and si mulated di stribution functions of 
Dehnenl d2OO0h andlFuxl d200lh . 



locity dispersion is cr u = 0.1 8vo- This plot is 
indeed in ve ry good agreement with p revious 
test-particl e (foehnenll200ri lFuxll200ll) and N- 
body (Fux 2001) simulations. 

By a quantitative comparison of the 
observed with the simulated distributions, 
iDehnenl (120001) deduced the Milky Way bar 
pattern speed to be Qb /^o = 1.85+0.15. In this 
calculation only the local velocity field, i.e., the 
observed u - v distribution, was taken into ac- 
count. Now, if one could also relate the dynam- 
ical effect of the bar to the derivatives of the ve- 
locities, an additional constraint on bar param- 
eters would be provided. Since the gradients of 
u and v are hard to measure directly, one needs 
an indirect way of achieving this task. The ob- 
vious candidates, as anticipated, are the Oort 
constants C and K. Oort's K, however, is hard 
to measure as mentioned above; hence only C, 
as estimated by OD03, will be employed here. 
By Fourier expansion of Eqs.[T](see Paper I for 
details) we estimate C from our numerically 
simulated velocity distributions for the initially 
cold and hot disks. 



4.1. Variation ofC with bar pattern 
speed and orientation 

In Fig.|2]we present our results for C as a func- 
tion of the bar angle, 0o (the angle by which the 
Sun's azimuth lags the bar's major axis). Each 
column shows a simulation with a different pat- 
tern speed, indicated in each plot. Rows from 
top to bottom show C as calculated from sam- 
ples at average heliocentric distances corre- 
sponding to d — 200, d = 400, and d = 600 pc, 
for a Solar radius rrj = 7.8 kpc. Solid and dot- 
ted lines represent the results for cold and hot 
disks, respectively. The dashed lines indicates 
C — 0. C is presented in units of Qo = vqItq. 
To make the discussion less cumbersome, we 
write Ch and C c to refer to the values for C as 
estimated from the hot and cold disks, respec- 
tively. 

Ch (dotted lines in Fig.|2j varies with galac- 
tic azimuth as Ch(0o) ~ sin(20o) for all of 
the fib values considered. On the other hand, 
the cold disk values (solid line) exhibit dif- 
ferent variations, depending on Q.^ or equiva- 
lently, on the ratio to/Volr- Closer to the OLR 
(left columns of Fig. |2j, C c (<po) approaches 
the functional behavior of Ch(0o)- Away from 
the OLR C c ((po) is shifted by 90° compared to 
Ch(0o)- While both cold and hot disks yield in 
increase in the amplitude of C(0o) as the pat- 
tern speed nears the OLR, the effect on the 
C c {<po) is much stronger. This is consistent with 
our expectation that the cold disk is affected 
more by the bar, especially near the OLR. 
While close to the OLR \C h (<p )\ < |C c (0 o )l, 
we observe the opposite behavior away from 
it. This could be exp l ained by the results of 
Muhl bauer & Dehnenl d2003l) . where they find 
that high velocity dispersion stars tend to shift 
the 'effective resonance' radially outwards. 

According to OD03's 'mode-mixing' cor- 
rected value, C ~ for the cold sample and 
decreases to about -10 km s~' for the hot pop- 
ulation. Hence we look for locations in Fig. [2] 
complying with this requirement. In the transi- 
tion region, 1.83 < fib/fio ^ 1.91, where the 
function C c (0o) transforms into its negative, 
there are specific angles which provide good 
matches for the observations, i.e., C c (</>o) ~ 
while Ch(0o) is significantly negative. We 
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Fig. 2. Each panel shows the variation of the Oort constant C with bar angle (p , for a simulation with 
the parameters given in Table [T] and a particular bar pattern speed, Q.\,, and a mean sample depth, d. Solid 
and dotted lines correspond to cold- and hot-disk values, respectively. Columns from left to right, show an 
increasing Sit, in units of Oo. Note that the OLR is at Qolr ~ 1.7. Different rows present results from samples 
with different mean heliocentric distance, d. Good matches to the observed trend in C (vanishing value for 
the cold disk and a large negative for the hot one) are achieved for 20° < 4>o < 45° and 1.83 < flb/^o < 1.91. 



performed simulations with different pattern 
speeds in the range 1.7 < Clb/Qo < 2.5. It 
was found that only in the transition region in 
Fig.|2]could one achieve a satisfactory match to 
the observations. Although the left- and right- 
most columns are nearly consistent with our 
requirement on C, we reject these values of 
Qb for the following reason. OD03 estimated 
a large negative C from both the reddest main- 
sequence stars, which are nearby, and the more 
distant red giants. This is inconsistent with the 
largest mean sample depth shown in Fig.[2]with 
Q b = 1.81, 1.93. We conclude that the bar pat- 
tern speed must lie in the range 



Ob/Oo = 1.87 + 0.04. 



(5) 



The bar angle with respect to the Sun has 
been proposed, as derived from IR photome- 
try, to lie in the range 15° - 45° with the Sun 
laggin g the bar major axis, whereas iDehnenl 
(2000) found that the bar reproduced the u- 
anomaly for 10° < <po < 70°. Examining 
Fig. [2] we find that our requirements on C re- 
sult in constraining the bar angle in the range 
20° < 0o < 45°, depending on Q. h . This result 
for the dependence of pattern speed on bar an- 
gle is in a very good agreement with Fig. 10 by 
IDehnenl (l2000h which shows a linear increase 



of the derived Qb/^o as a function of <pQ. This 
is remarkable since Dehnen (2000) obtained 
his results in a completely different way than 
our method here. 

The Oort constants A and B were found to 
also be affected by the bar, although to a much 
smaller extent. After correcting for asymmet- 
ric drift, we found that for the cold disk the bar 
causes AA = 3.5 + 0.7 and AB = 3.0 + 0.6 
km s _1 kpc , while the hot disk yielded AA = 
0.5 + 1.5 and AB = ± 1 km s" 1 kpc -1 . These 
values are calculated for ro = 7.8 kpc and 
vo = 220 km s _1 . Despite the bar's effect on 
C, the hot population does provide better mea- 
surements for A and B. 



5. Conclusion 

We have shown that the Galactic bar can ac- 
count for a trend seen in measurements of 
Oort's C, namely a more negative C-value with 
increasing velocity dispersion. At the same 
time we have improved the measurement of the 
bar pattern speed, finding Qb/^o = 1.87+0.04. 
In addition, the bar angle is found to be in 
the range 20° < (po < 45°. O ur result for Q b 
lies well within the estimate b vlDehnenld2000l) 

on 



lies well witmn tne estimate by Dennen UU 
and that by bebattista et al.1 (120021) . based 
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OH/IR star kinematics. This study provides an 
improvement on the measurement of the bar 
pattern speed b y a factor of ~ 4 compared to 
previous work dDehnenl|2000l) . The improved 
constraints on bar parameters should be tested 
by, and will benefit, future studies of the bar 
structure in the Galactic center region, that 
will become possible with future radial veloc- 
ity and proper m otion studies (e.g., BRAVA, 
iRich et all2007l) . 

In addition to a flat RC we have also con- 
sidered a power law initial tangential velocity 
V <P = v o{ r l r of with/? = 0.1,-0.1 correspond- 
ing to a rising and a declining RC, respectively. 
We found that for j3 — 0. 1 no additional error 
in Q/, is introduced. However, for /3 = -0. 1 we 
estimated Q b /Qo = 1-85 ± 0.06. From Fig. [2] 
we see that the main source of error in the bar 
pattern speed is due to the u ncertainty in the 
bar an gle. A recent work bv iRattenburv et al] 
(120071) used OGLE-II microlensing observa- 
tions of red clump giants in the Galactic bulge, 
to estimate a bar angle of 24° - 27°. Using this 
as an additional constraint we obtain flb/fio = 
1.84 + 0.01. 

While we account for the trend of increas- 
ingly negative C with increasing velocity dis- 
persion, we fail to reproduce the size of the C- 
value measured by OD03. The most negative 
value we predict is about -6 while they mea- 
sure -10 km s _I kpc -1 . We discuss possible 
reasons for this discrepancy: (i) we expect that 
the magnitude of C is related to the fraction 
of SN stars composing the Hercules stream. 
It is possible that during its formation the bar 
changed its pattern speed forcing more stars to 
be trapped in the 2: 1 OLR, thus increasing this 
fraction. This would give rise to a more neg- 
ative C from hot stars, while leaving the cold 
population unaffected; (ii) the numerical model 
here uses a distance limited sample. Modeling 
a stellar population with a magnitude limited 
sample is more appropriate to compare to the 
observed measurements, (iii) OD03's 'mode- 
mixing' correction is only valid if the Sun's 
motion in the z-direction were assumed con- 
stant. This is invalidated if, for example, a local 
Galactic warp were present. 
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